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Abstract 

A discrete nonlinear system is analysed in case of open chain boundary conditions at the ends. 
It is shown that the integrability of the system remains intact, by obtaining a modified set of 
Lax equations which automatically take care of the boundary conditions. The same Lax pair 
also conforms to the conditions stipulated by Sklyanin [5]. The quantum inverse problem 
is set up and the diagonalisation is carried out by the method of sparation of variables. 
Backlund transformations are then derived under the modified boundary conditions using 
the classical r-matrix . Finally by quantising the Backlund transformation it is possible to 
identify the relation satisfied by the eigenvalue of Baxter's Q-operator even for the quasi 
periodic situation. 



1 Introduction 



Analysis of nonlinear integrable systems is a subject of immense importance both from 
the physical and mathematical points of view. The inverse problem formulated on the 
basis of Lax operators can solve a wide class of problems. On the other hand, it is also 
important that problems be formulated with pre-assigned boundary conditions and solutions 
of the corresponding Cauchy problem be obtained. However at times the imposition of 
finite boundary conditions may lead to a loss of integrability. An important developement 
in this regard was the discovery by Sklyanin [1], who showed how non-periodic boundary 
conditions could be imposed on an integrable model without destroying its integrability. In 
this communication we have analysed a discrete nonlinear system, (also known as the DST 
model), initially studied by Ragnisco et al and also recently by Sklyanin et al from the view 
point of its integrability in presence of open boundary conditions. Our approach is different 
from that of Sklyanin in the sense that we have changed the form of the Lax operator at 
the end points where the boundary conditions are imposed. We have then shown that our 
approach can also accomodate the conditions laid down by Sklyanin for the system to be 
integrable under such circumstances. The present formalism was first used by Zhou [3] for 
a fermionic spin system. Both the classical and quantum hamiltonians have been derived 
explicitly. The quantum R matrix has also been deduced so as to formulate the quantum 
inverse problem for the model under finite boundary conditions. Diagonalisation of the 
hamiltonian is carried out by using the method of separation of variables, as the boundary 
matrices are non diagonal in character. Of late there has been a great deal of interest in 
the study of canonical Backlund transformations within the framework of classical r matrix 
theory, specially under periodic boundary conditions [15]. We have discussed this issue in 
case of quasi -periodic and finite boundary conditions. Finally in the case of the less stringent 
quasi-periodic boundary conditions we have shown how a connection with the eigenvalue of 
Baxter's Q operator may be established [11] by quantising the Backlund transformation. 
The model under analysis is described by the following equations of motion. 

• _ 2 
Qn — Qn+l — Q n r n 



r n = -r„_i + q n r 2 n 
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where q n stands for the time derivative of q n . It was originally proposed by Ragnisco et al [2] 
in their analysis of Lie-Backlund symmetries of discrete systems. The Lax pair associated 
with (1.1) may be written as 



n+l 



L \& 



where 
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L n (\) = I i M B (A) 

r n -i y 

Note that consistency of (1.2) yields the equations of motion only for periodic boundary 
conditions. If however certain nontrivial boundary conditions are to be introduced then one 
has to adopt a different strategy as will be explained in the sequel. The system given by 
(1.1) possesses a hamiltonian structure with the following symplectic form [4]. 
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For the purpose of introducing the boundary conditions we modify the Lax operator in (1.2) 
in the following manner, 

= LjiX^j, j = 1, 2, ....TV 



1 -M j {\)^ j j = 2, 3, ....N — 1 
= W N+1 (\)V N+1 



dt '^'v 



dt 



~df 



(1.5) 



where W^+i, W\ are two new 2x2 matrices depending on the spectral parameter A and on 
the dynamical variables. The usual consistency condition of (1.2) viz 

dLj{\) 



dt 



M j+1 L,(A) - Iy(A)M,-(A) 



1.6) 



is now replaced by the following set of equations: 
dLj{\) 



dt 



M j+1 (\)Lj(\) - H\)M 3 (\) (j = 1, 2....JV - 1) 



dL N {\) 
dt 



= W N+1 L N (X) - L N (X)M N (X) 



w N+1 = 




A 

"2 J 
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= M 2 (A)Lx(A) - L 1 (X)W 1 (X) (1.7) 

That is the consistency conditions at the two ends are different. In the present case we have 
formally. 

( 4 0^, \ / 4 0l \ 

(1.8) 

so that upon imposing the boundary conditions r = 9-,qw + i = 6+ the equations for the 
two ends of the discrete chain turn out to be as follows. 

<ii = q2-qfr 1 fi = -(6L - q x rl) 

Qn = (0+ - q 2 N r N ) r N = -(r^-i - 5jv^) (1-9) 
We shall refer to the above system as an open chain. 

2 Boundary conditions and classical r— matrix 

To ascertain the hamiltonian associated with the open chain system it is always useful to 
compute the so called classical r-matrix through the poisson bracket given by the symplectic 
form in (1.4) [8]. It is straight forward to show that. 

{L n (X) ® ; L m {n)} = [r(A, /x), L n (X) <g> L n (fi)]5 nm (2.1) 

where r(A,/i) = —jz-^,'P standing for the permutation operator. The mondromy matrix is 
then defined by 

T N (X) = J] L n (X) (2.2) 

n=l 

From (2.1) and (2.2) it follows that 

{T N (X) ® : T N (fi)} = [r(X, fi),T N (X) <8) T N (/j,)\ (2.3) 

We shall refer to this poisson algebra as the classical inverse scattering method (CISM-I) 
algebra. It follows that 



dT, 



N 



N 



(A) = ]T L N ...L n+1 L n L n . 1 ....L 1 (X) = W N+1 T N (X) - T N (X)W 1 (2.4) 



from which it is easy to conclude that J(A) = trTTv(A) is a constant of motion in the periodic 
case i.e when q m+N = q m and r m+N = r m . 



\ 



c(0 



(2.8) 



2.1 Quasiperiodic case 

Instead of the periodic boundary conditions that are customarily assumed, one can impose 
the so called quasiperiodic boundary conditions on the system, without destroying its inte- 
grability. This requires the existence of a matrix C(A,£) such that 

J(\) = tr[C(\,t)T N (\)] (2.5) 

is time independent so that in the present situation one has 

= tr[C(\,0(M N+1 T N -T N M 1 )} = tr[C(\,£)M N+1 T N ] - tr^C^A, 0^] (2.6) 

from which it follows that 

C(\,OM N+1 = M 1 C(\,0 (2.7) 
As a specific example one can take 

o d J 

To analyse the structure of J(A) in detail, we note the following asymptotics of the mon- 
odromy matrix in our case. 

T 1 N 1 (\)^\ N + \ N - 1 S + p 2 \ N - 2 + .... 
T£(\)*\ N - 1 q 1 + 0(\ N - 2 ) 
T 2 N 1 (X)^X N - 1 r N + 0(X N - 2 ) 
T 22 ^\ N - 2 p' 2 + 0(\ N - 2 ) (2.9) 

whence we get 

J(A) = tr(C(OT N (X)) = r*A" + (r^S)X N - 1 + (r^A^ 2 + (^)A"" 2 + .... (2.10) 
where 

N N-l 

S = J2 S ii V2 = J2(<li+l r i) +J2 S i S j' P2= r Nqi (2.11) 
i=l 1=1 j<jr 

One can verify, using the equations of motion, that each coefficient of A J is a constant. 
The following nonlinear combination of these constants then gives the Hamiltonian in the 
quasiperiodic case. 

1 JV_1 1 

h = -»ch& - 2 P2 ) + eV 2 = r*[E fc+ifi - 2 E ^ 2 + ^ft] ( 2 - 12 ) 

1=1 I 



2.2 Classical open chain 

As mentioned earlier, the open chain in which independent boundary conditions are imposed 
at the two ends, represents a more general case compared to the quasi-periodic boundary 
conditions. The equations of motion are given by ( 1.1) and ( 1.9). Now by following the 
arguments given by Sklyanin in [1] one can show that 

r(A) = tr[K + (X)U(X)}, U(X) = T N (X)K ( X)T N \-X) (2.13) 

represents the generator of the infinite conserved quantities with the matrices K±(\) in- 
corporating the effects of the nontrivial boundary conditions. In our case the classical r— 
matrix being symmetric, both K± satisfy the same equation namely 

[r(A, //), K(X) <g> K(n)] + K\X)r(X + fi)K 2 (fi) - K 2 (fi)r(X + fi)K\fi) = (2.14) 

where K 1 (X) = K(\)®I and K 2 (X) = I®K(X). The following solutions of (2.14) are found 
to be admissible . 

I r7_i_ si 

K + (X) 



9_ X \ 



d_ 



J 



0+ 

a e A 



(2.15) 
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By explicit multiplication we find using (2.9) and ( 2.15) in ( 2.13) 

r(A) = X N+2 + 2X N (a 2 - ]-S 2 + 9l 0_ + r N 6 + ) + 0(A^ 2 ) (2.16) 

where a 2 = Ekj=i s i s j + Y.fji 1 Qi=i r 'i and S = X^Li s i so that the hamiltonian is given by 
the expression. 

N-l i N 

H=Y, Qi+m - -z £ <lM + <h°- + (2.17) 
i=i / i=i 

To ascertain if f (A) = we observe that 

dr 

-(A) = tr[K + (A)T JV (A)ir_(A)T w 1 (-A) + K + (A)T W (A)K_(A)(-T W 1 (-A)T W (-A)T JV 1 (-A))] 
dr 

— (A) =tr[^+(A)W^(A)C/(A)] - tr[K + (A)T JV (A)^ 1 (A)^(A)T/(-A)] 
at 

-tr[K + (A)C/(A)^(-A)] + tr[K + (A)T w (A)K_(A)Vy 1 (-A)T w 1 (-A)] (2.18) 
hence f(A) = if and only if 

K + (X)W N+1 (X) = W N+1 (-X)K + (X) (2.19) 



Wi(\)K-(\) = tf_(A)Wi(-A) 
which are found to be identical to the conditions deduced by Sklyanin [5]. 



(2.20) 



3 Backlund transformation 

In this section we derive a Backlund transformation (BT) for the model under consideration, 
using the recently developed approach of Sklyanin [6, 7]. This relies on the algebra of the 
local lax operator as given by (2.1). For the sake of clarity we shall begin by retracing the 
methods used by Sklyanin in [16] and then consider the case where the nonlinear variables do 
not obey periodic boundary conditions. Sklyanin's formalism requires us to look for another 
representation of the local algebra given in (2.1). Let g^X — a) be such a representation 
satisfying 

- a) <g>, gi{n - cr)} = [r(A - //), #(A - a) <g> gfa - a)] (3.1) 

where A and \i are spectral parameters and cr is a complex parameter of the transformation. 
Moreover gi{X — cr) is assumed to depend on a canonical set of variables (ij,Tj). In (3.1) 
r(A — ji) = —-^ L - To construct the BT 

B a : T N (x, X- A) — T N (y, Y; A) (3.2) 

we shall first construct a local transformation 

bf -.L^Xi-X) -^L^Yi-X) (3.3) 

defined by the following 

gi(X - cr; S;, S , i )L i (A; x h X{) = L;(A; y h Y^g^X - a; U, %) (3.4) 

here Li(X;Xi,Xi) is the same as in (1.3), we have simply relabelled the variables qi,r{ by 
Xi,Xi for notational convenience. In (3.4) we see that in addition to the canonical variables 
(xi,Xi)f =1 and (yi,Yi)f =1 defining the BT we have two additional sets of canonical variables 
(siiSi)^ and (U,Ti)f =l . It is assumed that all the four sets of canonical variables are 
mutually commuting quantities. The transformation b$ is therefore a transformation in an 
extended phase space. Now by imposing suitable constraints it is possible to eliminate the 
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Qi(X - a); Si, Si) = 



V 



(3.5) 



auxiliary variables (si,S — i) and (U,Ti) from (3.4) thereby obtaining a local transformation 
We get from (3.4) the following set of independent equations, taking ^(A — a; s i: Si) to 
be of the form. 

1 

~Si A — a — SiSi 
(xi + Si)Xi = Di(Yi - Ti), -SiXi(xi + =Yi-Ti + aX { yi 

U = -yi, Xi = Si (3.6) 

Imposition of the constraints Sj = U+i and Si = T i+ i in these equations leads to 

X i = -- — Yi = X^ + ^^X t (3.7) 

Vi Xi - y i+1 yi 

which defines the local transformation b$ : (xj, Xi) — > (yi, Yi). It should be pointed out that 

this local BT may be obtained from a generating function G a where 

N _ 

G CT = E [^—^ + o log(^^±i) ] (3.8) 

i=l Vi Vi 

X, = Y t = f^ (3.9) 

oxi dy t 

Moreover as a consequence of the constraints s$ = t i+ i, Si = T i+1 and (3.6) we see that 

gi(\ - a; s h Si) = g i+1 (X - a;U,Ti) (3.10) 

and hence (3.4) becomes 

g i+ i(\ - a;y i+ i, Xi)Li(X; Xi, Xi) — Lj(A; XiXi)gi(X - a; y iy X^i) (3.11) 

Now it is well known that the generating function for the integrals of motion are given by 
the transfer matrix J(A) defined as the trace of the mondromy matrix, i.e 

J(A) =trT N (x,X;\) (3.12) 

where T N (X) obeys the CISM-I algebra. 

{Tjv(A) <g> T N {n)} = [r{X-n),T N {X)®T N {ii)\ (3.13) 

Under the BT defined by (3.11) the transfer matrix assumes the form. 

J(A) = tr[g N l +l (X - a;y N+1 ,X N )T N (y,Y; A)^i(A - a;yi,X Q )} (3.14) 

7 



Clearly invariance of the integrals of motion and thus of the hamiltonian generated by J(A) 
requires that the variables in gi{X — a) be periodic i.e 

Xn+tyi = X m , UN+rn = Dm (3.15) 

In such cases the BT will be canonical in nature since the hamiltonian is invariant. 
In the quasi periodic case the boundary conditions on the dynamical variables are determined 
by (2.7), so that upon using the form of C(£) as given in (2.8) we find that the variables in 
Li (A) obey 

x N+1 =£X N , X = £X N (3.16) 

We refer to these as the quasi periodic boundary conditions. Now in view of (2.5) which 
defines the generator of the integrals of motion J(A) we see that as a consequence of the BT 
given by (3.12), it is modified to the following. 

J(A) = tr[C(X, 0^+i(A - a; y N+1 , X N )T N (y, Y; X) 9l (X - a; Vl , X )\ (3.17) 

Clearly invariance of J(A) requires that 

<7i (A - a; yi , X )C(X, 0^+i(A - a- y N+1 , X N ) = C(X, £) (3.18) 

which gives. 

y N+1 = £ yi , X = £X N (3.19) 

Structurally these equations are consistent with (3.16) because the Lax operator Lj(A; y^ Yi) 
in (3.11) is obtained simply by replacing the variables (x i: Xi)f =1 with (yi,Yj)f =1 . Therefore 
the BT generated by (3.11) through the matrix gj(X — a) is a canonical transformation 
even when the nonlinear dynamical variables obey quasiperiodic boundary conditions. It 
is now natural to enquire if a BT can be constructed in case of open chains. Infact the 
most general type of boundary conditions for Liouville integrable lattices was formulated 
some years back by Sklyanin in terms of the representations of a new algebra known as the 
reflection equation algebra (REA). The classical limit of the REA is given by the following 
Poisson algebra, which shall be referred to as the CISM-II algebra. 

{T^AXT^/i^KA-zi),^ 

(3.20) 
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Now if K±(\) be representations of the above algebra then it may be shown that U(X) = 
T N {X)K^{X)T N \-X) also satisfies the CISM-II algebra, where T N {\) is defined by (2.2). 
Consequently it follows that 

r(A) = trK + (X)U(X) (3.21) 

is a generating function of the integrals of motion. We have already seen that under the BT 
given by (3.11) the monodromy matrix undergoes the following transformation. 

B a : T N (x, X; A) — ► S^+^A - a; y N+1 , X N )T N (y, Y; X)gi(X - a; y u X ) (3.22) 

Next consider 2x2 matrices V±(A) and construct the following quantities: 

Z(X) = T N (x, X; X)V4X)T N \x, X; -A) J(A) = tr(V + (X)Z(X)) (3.23) 

note that Z(X) per se might not satisfy the CISM -II algebra. But if we apply a BT as given 
by (3.11) then J (A) assumes the following form. 

J(A) = tr[{g N+1 (-X - a)V + (X)g^ +1 (X - a)}T N (y, Y; A)x 

x{ 9l (X - a)V4X)g^(-X - a)}T^(y, Y; -A)] (3.24) 

It is evident that J(X) would be a generator of the integrals of motion provided that the 
quantities in curly brackets are representations of the CISM -II algebra. This infact opens up 
the possibility of deriving a hierarchy of suitable boundary conditions for integrable systems 
starting from appropriate matrices V + . We will examine this possibility by deriving the 
usual scalar boundary matrices. First of all if the quantities in curly brackets were similarity 
transformations, then we could have taken the matrices V± to be simply the matrices K±(X) 
as defined in ( 2.15). This is a consequence of the observation that the CISM-II algebra 
is invariant under a similarity transformation of its elements. However this possibility is 
spoilt by the fact that the arguments of gi( N+ i) and 5^+1) i n ( 3.24) are different. Consider 
therefore the following forms of the matrices V±(X) 

\ 



X I a+4 b 



V + (X) 



x 



X + \ 1 -a + { j 



(3.25) 
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X 



X-a 



XA 1 + A + { (3X 2 + B 



so that 



and 



g N+1 (-X - a)V + (\)gJ? +1 (\ - a) = K + (X) = 



g 1 (X-a)V-(X)g^(-X-a) = K + (X) 



fe + o ^ 
yx e +J 



9_ X 
9_ 



J 



Here 



a = y N+1 -9 + d=-a9 + , b = y N+1 (29 + - y N +i) 
and A 1 = X , A = -(aX + 9_ - Vl X 2 ), 5 = o9_ 

13=1, C = X 2 , B = -( yi X -a) 2 + 29_ yi 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



Disregarding the overall factor in V + (X) we see that its structure is similar to the general 
scalar solution of the CISM -II algebra. While that of V-(X) is similar to that derived in 
[10] ( compare their eqns ( 2.3) -(2.5) with a — 7 = 0, f3 — 1). The difference in our case 
is that all the entries in V_(A) are scalars. Thus while the former is by itself a solution of 
the CISM -II algebra, for V_(X) to do the same it is necessary that X = 9_ = 0. The other 
pertinent issue is whether the matrices V±(A) are compatible with the integrability of the 
model under consideration, the constraint of integrability requires these matrices to satisfy 
(2.20) which is a stringent condition. In case of V+(A) the boundary condition yjv+i = 9 + 
ensures a = b = so that the integrability condition is not fulfiled unless Y/v = constant 
which clearly is in contradiction to the fact that (y^, Yi)f =1 are dynamical variables. Thus 
Kt(A) by themselves cannot be made to satisfy the requirements of integrabilty and CISM 
-II algebra simultaneously. However as a result of the BT the quantities Z(X) and J(A) are 
modified in a such a manner that the new boundary matrices K±(X) appearing in lieu of 
V±(A) are not only solutions of CISM -II algebra but are also compatible with the integrability 
of the model under consideration. 
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4 Quantum open chain 

In this section we consider the quantum inverse problem for a open chain. Note that the 
Lax operator in (1.3) is essentially similar to the discrete self trapping (DST) model [16]. 
The integrability and quantum inverse problem in absence of finite boundary conditions were 
studied in [12, 13]. From the hamiltonian structure of our system we have seen that {q n -,r n ) 
are a conjugate pair of variables; so in the quantum version of the model we impose the 
following commutation relations. 

[q n , r m ] = V$nm (V = ity [in, q m ] = [r n , r m ] = (4.1) 

As before the monodromy matrix is defined by T N (X) = IT^Li L n (X) and obeys the QISM -I 
algebra [8] 

R 12 (X - M(X)T%(jj.) = 7* (A)fli 2 (A - ( 4 - 2 ) 

with Ru(X — fj) = 1 + ^r'P- The matrices determining the boundary conditions K±(\) 
obey the conditions [9]. 

R 12 (X - fi)K w (X)R l2 (X + fi)K (2) (X) = K^(i2)R 12 (X + ^)K { l ] R l2 (X - /i) (4.3) 

R 12 (-X + p)K^(X)R 12 (-X -ii- 2r ) )K (2)t2 (p) 

= K^ h (fj,)R 12 (-X -ii- 2r ) )K^ h (X)R 12 (-X + p) (4.4) 

The form of K±(X) given in (2.15) is seen to satisfy these equations. The generators of the 
conserved densities are then built from. 

/ 



U(X) 



A(X) B(X) 
C(X) D(X) 



T(X)K^(X — 77/2, ^)a 2 T(—X)a 2 (4.5) 



which obeys the algebra 

R 12 (X - fi)U l (X)R 12 (X + - v)U 2 {X) = U 2 {X)R l2 {X + /1 - r ] )U 1 (X)R 12 (X - fj.) (4.6) 
and is given by 

r(A) = tr[K + {X - rj/2, £ + )[/(A)] = X 2N + 2 + 2X 2N H q + (4.7) 
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where the hamiltonian with quantum corrections is given by 

JV-l ^ N J2 

2 



1 n 

H q = E ^+1^ - o £(^) 2 - -5- + + £-<?i (4-8) 
i=i 8 



On the other hand with U(X) given by (4.5), we can write the general form of r(A) as [10] 

r(A) = e+(A(A) + £>(A)) + (A + ij/2)fl(A) (4.9) 

4.1 Separation of variables 

From (4.6) we can find the operator algebra for the elements A(X),B(X) etc. It is however 
more convenient to define instead of D(X) the following operator. 

D*(X) = 2XD(X) - r]A(X) (4.10) 

The commutation relations relevant for our purposes are then obtainable from (4.3) and are 
as follows. 

[fl(A),fl(/i)]=0 (4.11) 

(4.12) 

Furthermore from the definition of £/(A) as given in (4.5) the operators A(X),B(X) etc are 
found to have the following forms: 

A(X) = (-) N r N {X 2N + (S + r 1 /2)X 2N ~ 1 + } (4.14) 

D(X) = (-) N r N {X 2N + (S- r 1 /2)X 2N - 1 + } (4.15) 

N 

5(A) = (-f(X - V /2){X 2N + } = (-)"(A - v /2) J] (A - u a )(X + u a ) (4.16) 

a=l 

Where S = Y^=iQi^i- From ( 4.9) we note that the presence of the operator B(X) in the 
expression for the transfer matrix r(A) prevents us from setting up the usual algebraic Bethe 
ansatz. Consequently we have recast B(X) in the form given in (4.16) where u a are the zeros 

12 



of this operator. Henceforth {u a \ a = 1,2...N} are themselves to be regarded as operators 
defining B(X) and it follows from their commutativity given by (4.11) that 

[u a ,M = Q ( 4 - 17 ) 

We next define operators by left substitution A u a into the operators A(A) and D*(X) 
as shown below : 

v+ = D*(X^u a ), v~=A(X^u a ) (4.18) 
It then follows from the algebra of the operators A(A), D*(X) that 

[v±,Up}=±r)V±6 a p [*t,tf]=0 K,^]=0 (4.19) 

Notice that another set of dynamical variables can be constructed from the zeros of B(X) at 
A = —u a , ot — 1, 2...N and obeying 

[w±,Uf,] = ±r]w^5 a(3 [w+,w^] = [w-,wp] = (4.20) 

where tu+ = A(A =>- —u a ) and w~ = D*(X =>- —u a ). Then by using Lagrange interpolation 
one can rewrite the operators A(A) and D*(X) in the following manner . 

2N 2N \ _ \ 2N 

a(a) = e n y-^f-Aj+H^n^-Aa) (4-21) 

a=l P=l,f3^a *a — a=l 
2N 2N \ _ \ 2N 

d*(x) = j: n ^f-^+H^n^-Aa) (4.22) 

a=l 0=l,p^a — Xp a =l 
X a = U a , 1 < a < N 

X a = -u a N + 1 < a < 2N 

i<«<iv 

= N + 1 < a < 2N 

The transfer matrix then assumes the form: 

2u a r(u a ) = i + (2u a + rj)v- + £+v+ (4.23) 
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Now by extensively using the preperties of quantum determinants [14, 9] it can be shown 
that from the commutation relations for the operators u a ,v^,w^ one can define a Hilbert 
space where the eigenvalue problem reduces to 

2u a r(u a )*(tti, u N ) = £+{2u a - 77) A~(w a )^( ,u a - 77, ....) + £+A + (w Q )^(..., u a + 77, ....) 

(4.24) 

where 



A-(u) = £_ + („- |) A + (u) = (2« - 77)(e- - (« + |)) (4.25) 

Setting 



AT 

*K..%) = II^a) (4-26) 

a=l 

we have finally 

2w a r(w Q )0(w Q ) = £+(2w a + ?7)A~(w a )0(w a - 77) + £+A + 0(w a + 77) (4.27) 



5 Relation with Baxter's Q-operator 

In this section we consider a quantum mechanical version of the Backlund transformation 
and show how in the quasi periodic case it may be connected to Baxter's Q-Operator. For 
this purpose we quantise the dynamical variables occurring in the Lax operator given in (1.3) 
which then asssumes the form. 

\ 

(5.1) 



Li(\) 



* ~ VQi9 qi qi 
-Vd qi 1 J 

where we have set r\ = —r)d qi and rj] = rjSij. The quantum analog of the BT for integrable 
DST model which essentially has the same Lax operator as ours was studied in great detail 
by Kuznetsov, et al in [15]; others have also studied the same model from alternative points 
of view [16, 11]. The essential idea is to reproduce an operator satisfying the following 
conditions, 

J(X)Q(X) = A+(A)Q(A + 77) + A_(A)Q(A - 77) (5.2) 
Q(A)Q(/i) = QMQ(A) (5.3) 
j(»)Q(\)=Q(\)j(») (5.4) 
14 



and having a common set of eigenfunctions 

Q(A)$ = Q(A)$, J(A)$ = J(A)$ (5.5) 

with the eigenvalues satisfying a finite difference equation. 

r(A)Q(A) = A+(A)Q(A + 77) + 5_(A)Q(A - 77) (5.6) 

Incidentally the last equation arose in Baxter's analysis of the XXX model where r(A) was 

the transfer matrix of the model with Q(X) being a suitable polynomial. 

To prove (5.3) it is sufficient to demand the existence of an operator R a such that 

G{\ - a; t, dt)L(\; y, d y )R a = R a L(\; q, d q )Q{\ - a; s, d s ) (5.7) 

where we have dropped the indices V on the variables for brevity. Note that here Q(X — a) 
is assumed to satisfy the following algebra. 

Ri2(X - n)Q\X)Q\n) =)G 2 (ri)G\\)Ri2(\ - ») (5.8) 

and as before 

R(X — n) =1 + — ^ — V (5.9) 
A — jj, 

regarding the explicit form of Q(\ — a) we take it to be as follows. 

Q(X - a; s,d s ) = 



A — a — sd* s ^ 



(5.10) 



-7]d s 1 j 

On the basis of the above remarks proceeding as in [15] we finally arrive at a form of Q{a) 
as an integral operator with a kernel given by 

N 

2( (J ) = II^^^+i,!/!,?;) (5-11) 

i=i 



where 



Wi(<T;yi+i,yi,q i ) = 7r r(<T/ri + l)y i ( ) i exp{ ) (5.12) 



To derive the last property of Baxter's Q-operator we note that J(A) = tr[C(\,£)T N (\)] in 
the quasi periodic case. So we consider 

[Q(a)j(a)$](y) = tr j dq N Q{a)T N {a)<b{q) (5.13) 
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which yields for the kernel of the joint operator Q(a)J(a) 



N 



N 



[Q(a)J(a)$](y \ x) = tr{C{£)T* N {<j) «>»} = tr{C(£) JJ L*(a) Wi } (5.14) 



i=i 



i=i 



where L*(a) represents the adjoint of the corresponding local Lax operator. The last equation 
may be recast into the following form 



with 



[Q(<j)j(<j)<S>](y\x) = Q(a)tr[C(Of N (a)} 



r tt f / \ t~t ( <T + V + VQid qi log Wi qi 

Tn(<7) = 11 Li{&) = [[ | 

r)d qi \ogw t 



i=i i=i y r)O qi log Wi 1 

It is now possible to triangularise the matrix Lj(cr) by the following transformation: 

\ 



(5.15) 



(5.16) 



Li(cr) — > S , i+ 1 1 L i ((j)S'i, Si 



1 Vi+i 
1 



(5.17) 



This transformation leaves the quasiperiodic boundary matrix C(£) invariant, while causing 
Li(a) to become triangular as shown below. 



1^) 



It is then straightforward to show that 

TV 



(7- 





+1) 



(5.18) 



AT 



JV 



o n ^(-) = r ^ n - 1) + ^ n + 1) 

i=i ^ 



(5.19) 



i=l i=l 

This relation ought to be compared with the eigenvalue deduced by the standard method of 
algebraic Bethe ansatz for the quasi periodic case which is of the form: 



in in in 



(5.20) 



i=l 



i=l 



i=l 



6 Discussion 

We have analysed a discrete nonlinear integrable system in the presence of finite boundary 
conditions and have shown how the diagonalisation problem can be handled when the bound- 
ary matrices are upper (lower) triangular in nature. The construction of Backlund transfor- 
mations under different types of boundary conditions was analysed and for the quasiperiodic 
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case we have shown that a connectionbetween the quantised form of the BT and Baxter's Q 
operator may be established. 

One of us (AGC) wishes to acknowledge the financial support provided by UGC ( Govt, of 
India) PSW-026/99-00. 
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